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Abstract. 

Studies of the Interacting Boson Approximation (IBA) model for large boson numbers have 
been triggered by the discovery of shape/phase transitions between different limiting symmetries of 
the model. These transitions become sharper in the large boson number limit, revealing previously 
unnoticed regularities, which also survive to a large extent for finite boson numbers, corresponding 
to valence nucleon pairs in collective nuclei. It is shown that energies of 0+ states grow linearly with 
their ordinal number n in all three limiting symmetries of IBA [U(5), SU(3), and 0(6)] . Furthermore, 
it is proved that the narrow transition region separating the symmetry triangle of the IBA into 
S ' a spherical and a deformed region is described quite well by the degeneracies £(0^) = E{6'^), 

£■(0^) = E{10+), E{0^) = ^CMl), while the energy ratio £'(6|)/£'(0j) turns out to be a simple, 
empirical, easy-to-measure effective order parameter, distinguishing between first- and second-order 
transitions. The energies of 0+ states near the point of the first order shape/phase transition between 
U(5) and SU(3) are shown to grow as n(nH-3), in agreement with the rule dictated by the relevant 
critical point symmetries resulting in the framework of special solutions of the Bohr Hamiltonian. 
Xf-s [ The underlying partial dynamical symmetries and quasi-dynamical symmetries are also discussed. 
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The Interacting Boson Approximation (IBA) model [1], describing collective phe- 
nomena in atomic nuclei in terms of s and d bosons (of angular momentum and 2, 
respectively), is known to possess an overall U(6) symmetry, containing three different 
c^ I dynamical symmetries, U(5), SU(3), and 0(6), corresponding to near- spherical (vibra- 

tional), axially symmetric prolate deformed (rotational), and soft with respect to axial 
asymmetry (y-unstable) nuclei respectively. These limiting symmetries are shown at the 
vertices of the symmetry triangle [2] of the model, shown in Fig. 1 . 

An energy functional can be obtained M in the classical limit of the model, through 
the use of the coherent state formalism o, \&i Studying this energy functional in the 
framework of catastrophe theory one can see IM that a first order phase transition (in the 
Ehrenfest classification) is predicted to occur between the limiting symmetries U(5) and 
SU(3), while a second order phase transition is expected between U(5) and 0(6). We 
refer to these transitions as shape/phase transitions. A narrow shape coexistence region 
is then predicted [4] in the symmetry triangle of the IBA, separating the spherical and 
deformed phases. The shape coexistence region shrinks into the point of second order 
phase transition as the U(5)-0(6) line is approached, as shown in Fig. 1. 
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FIGURE 1. IBA symmetry triangle with the three dynamical symmetries. The critical point models 
E(5) and X(5) are placed close to the phase transition region (slanted lines). The solid curve indicates the 
AUiassid-Whelan arc of regularity. Adopted from Ref. ^. 



Shape/phase transitions have been considered recently in the framework of the geo- 
metric collective model [8] as well, in which the collective variables /3 and 7 are used. 
The second order transition between U(5) and 0(6) has been described by the E(5) crit- 
ical point symmetry [9], using in the Bohr Hamiltonian [10] a potential independent of 
7 and having the shape of an infinite square well potential in j8, while the first order 
transition between U(5) and SU(3) has been described by the X(5) critical point symme- 



try nun , using in the Bohr Hamiltonian a potential of the form m(/3) +v{y), with m(/3) 
having again the shape of an infinite square well potential in /3, and v{y) being a steep 
harmonic oscillator centered around 7 = 0. E(5) and X(5) are shown in Fig. 1, close to 
the points of the second and first order phase transitions of the IBA, respectively. 

The IBA model predictions have also been studied using different measures of chaotic 
behavior. It has been found that chaotic behavior prevails over most of the symmetry 
triangle, with the noticeable exception of two regions, in which highly regular behavior 
is observed [12]. One of them is located along the U(5)-0(6) line, and its existence is 
due to the 0(5) symmetry known to survive along the whole line [13]. The other regular 
region is connecting U(5) and SU(3) through a narrow path within the triangle, called the 
Alhassid-Whelan arc of regularity [12. ] (also shown in Fig. 1), the symmetry implying 
its existence being yet unknown. 

States with zero angular momentum are of particular interest, since centrifugal effects 
are absent, facilitating the detection of underlying symmetries. In the U(5) limit of 
the IBA, the energies of 0^ states increase linearly with the number of d bosons, 
corresponding to their phonon number. In the SU(3) limit, the energies of 0"'" bandheads 
are determined by the eigenvalues of the second order Casimir invariant of SU(3), the 
results shown in Table 1. In the 0(6) limit, they are determined by the second order 
Casimir invariant of 0(6), the results also shown in Table 1 . From Table 1 it is clear 
that in the limit of large boson numbers Nb, the energies of 0+ bandheads in all three 
dynamical symmetries of the IBA grow linearly, E =An. 

IBA numerical calculations have been performed using the recently developed code 



TABLE 1. Irreducible representations (irreps) of SU(3) and 0(6) and the corresponding 
energy of the excited 0^ bandheads. A^ stands for the boson number, Nb- Adopted from 
Ref.ll. 
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FIGURE 2. (Left) Line of degeneracy between the Oj and 6^ levels for A^b = 10, 40, 100, and 250 
in the IBA triangle. (Right) Line of degeneracy between the 0^ and lOj'^ levels for Nb = 250 (top) and 
between the 04^ and 14^^ levels for A^^ = 250 (bottom) in the IBA triangle. The dashed lines denote the 
critical region in the IBA obtained in the large Nb limit from the intrinsic state formalism. Adopted from 
Ref. klTJ. 



IBAR llJ, [iSD, which can handle large boson numbers. The standard two-parameter 
IBA Hamiltonian [ J16l ] has been used, depending on the parameters ^ and X- Examining 
the degeneracy £(0^) = E{6^), a hallmark of the X(5) critical point symmetry, we find 
rtlTI] that its locus in the IBA triangle is a straight line, falling withing the coexistence 
region for large boson numbers, as seen in Fig. 2. Similar results are obtained for the 
degeneracies E{0^) = E{10^) and E{0^) = E{14^), also shown in Fig. 2. 

Plotting the ratio £'(6]'^)/£'(0j) lIlTll . related to the first of the degeneracies mentioned 
above, we see in Fig. 3 that it exhibits the behavior expected [18] for an order parameter 
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FIGURE 3. The ratio £'(6i^)/£'(0^) as a function of C, for three values of x for (a) Nb = 15 and (b) Nb 
= 100. The inset to (a) shows the corresponding behavior for the order parameter V2 lUSll . Adopted from 
Ref. lEl. 
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FIGURE 4. (a) Experimental £'(6[^)/£'(0j) ratio as a function of neutron number for the Nd, Sm, Gd, 
and Dy isotopes, (b) Same for the Xe and Ba isotopes. For smaller neutron numbers, the Oj^ state was 
taken in the ratio if its B{E2) decay was consistent with the a = N — 2 state. This corresponds to A^ = 74 
in Xe and A^ = 76,78 in Ba. Valence (hole) neutron number increases to the left. Adopted from Ref. OITI] . 



of a first (second) order phase transition for ;t = —1.32 (x =0)- It is remarkable that 
experimental data around the N — 90 isotones, which are known to be the best empirical 
examples of X(5) [19i|20i|2l|,|22i|23T, shown in Fig. 4, exhibits a behavior very similar 
to the one expected for a first order phase transition, while data around ^-^^Ba, the best 
example of E(5) [|24l l23ll. shows the behavior expected for a second order transition. 

Now we turn attention to energies of 0+ states within the critical point symmetries 
E(5) and X(5), mentioned above, as well as within Z(5) [25], a similar solution of 
the Bohr Hamiltonian, also using an infinite square well potential in /3, with y^ Tz/6. 
As seen in Table 2, although the relevant energies look different in each model when 
normalized to the energy of the 2^^ state, they become identical if normalized to the 
energy of Oj [3]. Moreover, they increase as n{n + 3). Similarly, energies of 0+ states 
within Z(4) [26] [similar to Z(5), but with 7 fixed to Tz/6] increase as n{n + 2.5). Within 



TABLE 2. (Left) Energies of 0+ states in the E(5), Z(5), and X(5) models. Energies on 
the left are in units of £'(2['^) = LO, while in the column Norm, in units £(0^) = LO. The 
normalized results are identical for each of the models. The column IBA-Norm gives the 
normalized 0+ energies for a large Nb IBA calculation near the critical point. (Middle) Same 
for the Z(4) model. (Right) Same for the X(3) model. Adopted from Ref. H. 
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TABLE 3. Order v, dimension, D, of the model space and v for J'' = 0+ states in the geometrical 
models E(5), X(5), Z(5), Z(4), and X(3). J is the spin of the level, T — J/2, and riw is the wobbling 
quantum number |8] which is zero for 0+ states. 
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X(3) [|27|] [similar to X(5), but with 7 fixed to 0], they increase as n{n + 2), as also 
shown in Table 2. These results can be easily interpreted [3] by taking into account 
the order of the Bessel functions appearing as eigenfunctions in these models, given 
in Table 3, as well as the fact that the spectra of the Bessel functions Jy are found 
to increase approximately as n{n + v + 3/2), this result being exact only for v = 1/2 
^M- Consideration of the second order Casimir operator of E(n) [|28U26|] . the Euclidean 
algebra in n dimensions, shows [3] that the present situation is a partial dynamical 
symmetry f2^ of Type I [30], in which some of the states (the 0+ states in the present 
case) preserve all the relevant symmetry. The recent conjecture [31] of a partial SU(3) 
dynamical symmetry underlying the Alhassid-Whelan arc of regularity is also receiving 
attention. 

It is a nontrivial result [3] that the IBA near the critical point of the first order transition 
also yields energies of 0+ states increasing as n{n + 3), i.e. in the same way as critical 
point symmetries based on infinite square well potentials in 5 dimensions (degrees of 
freedom) predict, as shown in Table 2. 

In conclusion, 0+ states in the large boson number limit of the IBA exhibit many 
interesting properties, as well as degeneracies to non-zero angular momentum states, 
which invite further investigations into determining of the symmetries underlying these 
regularities. 
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